Improved chain mean-field theory for quasi-one-dimensional quantum magnets 
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A novel mean-field approximation for quasi-one-dimensional (QID) quantum magnets is formu- 
lated. Our new mean-field approach is based on the Bethe-type effective-field theory, where thermal 
and quantum fiuctuations between the nearest-neighbor chains as well as those in each chain are 
taken into account exactly. The self-consistent equation for the critical temperature contains the 
boundary-field magnetic susceptibilities of a multichain cluster, which can be evaluated accurately 
by some analytic or numerical methods, such as the powerful quantum Monte Carlo method. We 
show that the accuracy of the critical temperature of QID magnets as a function of the strength 
of interchain coupling is significantly improved, compared with the conventional chain mean-field 
theory. It is also demonstrated that our new approximation can predict nontrivial dependence of 
critical temperature on the sign (i.e., ferromagnetic or antiferromagnetic) of interchain coupling as 
well as on the impurity concentration in randomly diluted QID Heisenberg antiferromagnets. 



I. INTRODUCTION 

Space dimensionality plays an essential role in phase 
transitions and critical phenomena of quantum magnets. 
As the dimension is lowered, effects of thermal and quan- 
tum fluctuations generally become stronger. As a re- 
sult, the quantum antiferromagnetic Heisenberg model in 
two dimensions, for example, does not exhibit long-range 
order any more except at the ground state, ^ though a 
finite-temperature phase transition occurs in its three- 
dimensional counterpart. Furthermore, there is no long- 
range order even at zero temperature in one dimension.^ 
Such one-dimensional magnets arose great interests as 
many novel phenomena characteristic to systems with 
strong quantum fluctuations, e.g., Tomonaga-Luttinger 
liquid state or Haldane gap state, have been observed 
theoretically as well as in the experiments.— 

Real materials, however, can not be purely one dimen- 
sional but three dimensional, i.e., there always exist inter- 
actions between one-dimensional chains [interchain inter- 
actions) albeit it is much weaker, often by orders of mag- 
nitude, than the dominant interactions along the chains 
(intrachain interactions). Three-dimensional materials 
with strong spatial anisotropy are often referred to as 
quasi-one-dimensional (QID) systems. Indeed, in many 
QID materials a long-range order emerges at low tem- 
peratures, which is a direct consequence of three dimen- 
sionality of the system. In order to explain low-energy 
behavior of such QID materials correctly, a theory which 
properly incorporates the effect of interchain interaction 
is essential. 

So far, effects of weak interchain interactions in QID 
quantum magnets have been studied mainly by means 
of the chain mean-field approximation, where interchain 
spin fluctuations are ignored completelyi^i^ Within the 
framework of the chain mean-field approximation, a QID 
magnet is reduced to a single chain in an effective exter- 
nal field. For the latter system, fortunately there exist 
a couple of exact solutions, otherwise one can still use 
powerful analytic methods, such as bosonization, as well 



as numerical simulations, such as exact diagonalization 
or density-matrix renormalization group method, which 
are effective especially in one dimension. 

Naively, one may expect that the chain mean-field ap- 
proximation becomes more and more accurate not only 
qualitatively but also quantitatively, as the interchain in- 
teractions become weak enough compared to those along 
the chains. The recent theoretical studyS as well as the 
sensitive Monte Carlo simulations2i^ on QID spin mod- 
els, however, have revealed that this is not the case; there 
remains systematic error of the chain mean-field theory 
even in the weak interchain coupling limit. Instead, the 
critical temperature as a function of the interchain cou- 
pling is well described by a chain mean-field formula with 
a renormalized effective coordination number (or effective 
interchain coupling^^). Especially, in the QID classical 
Ising models, it is proved that the renormalization fac- 
tor is exactly given by the critical transverse field of the 
quantum phase transition of a two-dimensional quantum 
Ising model4 This analytic result demonstrates clearly 
that the weak interchain coupling limit of the QID mag- 
net is not the weak coupling limit, but is still the strongly 
correlated regime. 

The renormalization of effective interplane coupling is 
also observed in weakly-coupled two-dimensional planes 
For this quasi-two-dimensional (Q2D) system, a scaling 
theory was developedf^ in which it is expected different 
scaling behavior depending on whether the purely two- 
dimensional system is in a quantum critical regime or in 
a renormalized classical regime. This prediction has been 
verified by a recent quantum Monte Carlo simulation.J^ 

Thus it has become evident by the recent theoretical 
and numerical studies that we need a theory beyond the 
conventional chain mean-field approximation to describe 
the critical temperature of the QID magnets in the weak 
interchain coupling regime more accurately. Even worse, 
when the chain mean-field theory is applied to a sys- 
tem with quenched disorder, the random average in each 
chain is taken before the thermal average on the whole 
lattice. The adverse impact of interchanging the aver- 
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age operations is further nontrivial. It is naturally ex- 
pected that some of these drawbacks of the chain mean- 
field theory might be remedied by considering multiple 
chains, instead of a single chain, and taking the inter- 
actions between those chains into account in a proper 
way. Along this line, Sandvik has proposed a multichain 
mean-field theory and applied it successfully to the prob- 
lem of ground state staggered magnetization in a two- 
dimensional Heisenberg antiferromagnetiii 



In the present paper, we propose a different type of 
chain mean-field theory, i.e., the chain Bethe approxima- 
tion. Actually, it is well known that there are two dif- 
ferent formulations in traditional mean-field theories for 
classical spin models, that is, the Weiss approximatioi*i^ 
and the Bethe approximation. ^^'^"^ In the former, the ef- 
fective field is identified explicitly with the order param- 
eter. In the latter, on the other hand, the effective field 
is determined implicitly so that the local order parame- 
ter at the central spin and that of spins on the cluster 
boundary coincide with each other. In the Bethe mean- 
field theory, therefore, spin fluctuations between nearest- 
neighboring sites are taken into account even in the low- 
est order approximation. Furthermore, as the cluster size 
increases, the critical temperature by the Bethe-type ap- 
proximation converges to the exact value more rapidly, 
being free from logarithmic corrections which is observed 
in the Weiss- type theory. By applying the idea of the 
Bethe-type effective-field theory, we introduce an new 
chain mean-field theory for QID quantum magnets. 



The present paper is organized as follows: In Sec. 2, 
after a brief review of the conventional chain mean-field 
approximation, we formulate our Bethe-type chain mean- 
field theory, the chain Bethe approximation. In Sec. 3, 
the chain Bethe approximation is applied to the QID 
Heisenberg antiferromagnets, where we show that by us- 
ing our new approximation, the accuracy of the criti- 
cal temperature is significantly improved, compared with 
the conventional chain mean-field theory. In addition, 
we demonstrate that the new theory can predict a lower 
critical temperature for the ferromagnetic interchain cou- 
pling than in the antiferromagnetic case, which is also 
confirmed by the QMC simulation. In Sec. 4, we apply 
the chain Bethe theory to a random quantum magnet, 
the site-diluted QID Heisenberg antiferromagnet, where 
the existence of a finite critical impurity concentration, 
above which the long-range order does not emerge even 
in the zero-temperature limit, is predicted by using the 
chain Bethe approximation. The final section is for a 
summary and discussion. 



II. CHAIN BETHE MEAN-FIELD THEORY 

A. Conventional chain mean-field approximation 

The Hamiltonian of spin- 1/2 Heisenberg model on a 
QID simple cubic lattice is defined by 



i,j,k i,j^k 
i,3,k 



id,k 



(1) 



where S'ij-fc = (Sfj^k^ ^Ij.k^ S!,j,k) is an S = 1/2 quan- 
tum spin operator at site {i,j,k). We take the lattice 
2;-axis as the chain direction. The first term in Eq. ([T]) 
represents the intrachain interactions (with coupling con- 
stant J), while the second and third terms denote the 
interchain interactions {J'). In the present paper we con- 
sider only the case where the intrachain coupling J is an- 
tiferromagnetic (J > 0), whereas the interchain coupling 
J' is either antiferromagnetic (J' > 0) or ferromagnetic 
(J' < 0). Generalization to the case with ferromagnetic 
intrachain coupling is also straightforward. The last term 
in Eq. JT]) represents an external magnetic field conjugate 
to the order parameter, where the phase factor (j) is de- 
fined as (p — — sign(J'). 

In the conventional chain mean-field approximation, 
referred to as the chain Weiss theory hereafter, spin fiuc- 
tuations between the chains are ignored and replaced by 
an effective field, i.e., 



+ J'{Si^j,k) ■ S^'j'.k + const.. 



(2) 



where = {i±l,j) or {i,j±l). As a result, the orig- 

inal Hamiltonian ((!]) is decoupled into a set of indepen- 
dent chains [Fig. [Ha)] . The effective chain Hamiltonian 
for («, j) — (0,0) is written as 



He 



with Sk = S'o,o,fc and 
Mfc = (Si,o,/c) + (S'-i.o,fe) + (S'o,i,fc) + (S^o,-i,fc) 



(3) 



(4) 



In the low-temperature ordered phase, a finite mag- 
netization appears spontaneously even at /i = 0. We 
assume the magnetization is along the z-direction in the 
spin space, i.e., {Si-j^k) = (0, 0, (-l)'=TO(r)). The 
chain Hamiltonian ([3]) with /i = is then reduced to 

n. = jY,Sk- Sk+i - 4 I J'l m{T) Y,{~^fSl (5) 



Since the magnitude of the spontaneous magnetization 
does not depend on the position of the spin, the self- 
consistent condition 



m{T) = TOc(r,4|J'|TO(r)) 



(6) 
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(a) 



(b) 



FIG. 1: Clusters for (a) chain Weiss and (b) chain Bethe 
approximations. The solid and dashed lines denote the intra- 
chain (J) and interchain (J') interactions, respectively. In the 
chain Weiss approximation, an alternating effective field with 
the same amplitude applies on all spins. In the chain Bethe 
approximation, on the other hand, a staggered effective field 
applies only on the spins on the side chains (filled symbols). 



must be fulfilled, where 

mc(T,/.)^l^(-lf(5,-)e 



(7) 



is the staggered magnetization density of genuinely one- 
dimensional antiferromagnetic chain of length L. The 
average {• ■ ■)c in Eq. ([7]) means the expectation value of 
the one-dimensional chain (i.e. Eq. ^ with J' — 0), 
while (• • • ) [e.g. in Eqs. ^ and @] denotes the average 
with respect to the original (full-3D) Hamiltonian 

At high temperatures, on the other hand, no spon- 
taneous magnetization appears. Under the presence of 
small external magnetic field, however, a finite magneti- 
zation (S'jj-fc) = {Q,Q,(jf+^{-lfm{T,h)) is induced. In 
this case, the effective chain Hamiltonian ([5]) is modified 
as 



k+l 



[h + 4\j'\m{T,h)]J2{~l)'s:. 

k 

(8) 



Since m{T, /i) ^ 1 for /i <IC 1 in the disordered phase, 
one can consider only the lowest order in h: 



m{T, h)^[h + 4 I J'l m(r, h)]xciT), 



(9) 



where Xc{T) is the zero- field staggered susceptibility of 
one-dimensional antiferromagnetic chain. Noticing that 
m(T, h) in the both sides can also be written as hx{T) 
for /i <C 1, we obtain the following mean- field expression 
for the susceptibility: 



X{T) 



Xc(T) 



i-4|j'|xc(r)- 



(10) 



In terms of the chain Weiss theory, the critical temper- 
ature is thus given by the pole of the r.h.s. of Eq. (fTO| . 
For generic QID lattices the self-consistent equation is 
written as follows: 



l-zIJ'IXc(Te) =0, 



(11) 



where z is the coordination number of the lattice, i.e., 
the number of nearest-neighbor chains. Note that the 
self-consistent equation depends not on the sign of J', 
but only on its absolute magnitude. In other words, the 
conventional chain Weiss theory does not distinguish be- 
tween the antiferromagnetic and ferromagnetic interchain 
interactions. This is one of the major drawbacks of the 
conventional chain mean-field theory. Another problem 
is that most of physical quantities, such as the energy, 
specific heat, correlation functions, etc, are the same as 
those of the genuine one-dimensional chain at tempera- 
tures higher than Tc, since the effective field is propor- 
tional to the order parameter, which is zero at T > T^. 
We will see below that these disadvantages of the chain 
Weiss approximation are solved in the chain Bethe mean- 
field theory. 



B. Bethe- type mean- field theory 

In the chain Weiss approximation, the effect of inter- 
chain interaction was replaced by an effective field. The 
effective interchain interaction is directly related with 
the order parameter, and interchain spin fluctuations are 
thus ignored completely. The approximation can be im- 
proved by taking interchain spin fluctuations into account 
systematically. In the chain Bethe approximation intro- 
duced in this section, the interaction between the nearest- 
neighbor chains are taken into account exactly, and those 
around the multichain cluster are treated as an effective 
fleld. In the Bethe approximation ; ^'^i the effective field 
is determined so that the magnetization of the central 
spin and that on the cluster boundary coincide with each 
other. On the contrary to the Weiss theory, the order pa- 
rameter is not given explicitly by the effective fleld, but 
is its implicit function. 

Let us consider the chain Bethe approximation for the 
simple cubic lattice. In this case we prepare a cluster of 
flve chains [Fig. [T]^b)], where an effective fleld is applied 
only to the spins on the side chains (black circles). The 
effective Hamiltonian is written as 

4 4 
"^O = J ^ ^ Sa,k ■ Sa^k+l + J' ^ ^ S'o.fe ■ Sa^k 



a=0 fc 



a— 1 k 
4 



hY^{~lfSl, -{h + /.eff) ^ '^(-l)'^o,fc' 

(12) 



a=l k 



where a 
a = 1, • • 
(0,±1)]. 



denotes the central chain — (0, 0)] and 

,4 denotes the side chains = (±1,0) or 
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FIG. 2: ftoff-dependence of the local magnetization of the cen- 
ter spins (open symbols) and the boundary spins (filled sym- 
bols) of L = 64 chain Bethe cluster with J'/ J = 0.1. At 
temperatures higher than the critical temperature (squares) 
two curves intersect only at /leff = 0, while at lower tem- 
perature (diamonds) nontrivial solutions corresponding to the 
symmetry-broken phases appear. The dashed lines denote the 
tangent of each curve at h^s — 0. In the inset, temperature 
dependence of An [Eq. (|16|) ] is also presented. 



In the chain Bethe approximation, we impose the con- 
dition that the absolute value of the local magnetization 
does not depend on its position: 



for any fc, fc', and a — 1, - ■ ■ , 4, or equivalently 



(13) 



each other (Fig. [2])- In practice, we set ft, = in Eq. (fT4|) 
and expand the both sides in terms of hes- The self- 
consistent equation for the critical temperature is then 
written as 



a=l 



= 



(16) 



with the boundary-field susceptibilities 

dmfi^a{T,h,hcB) 



dhcs 



(17) 



7 X! (■^g.fci'S'a'.fcOf 



a' = l k,k' 



where {A;B)^-^ denotes the canonical correlation: 
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Tre 



(18) 



of two operators A and B. In the inset of Fig. [21 we 
show the temperature dependence of An (T) . As the tem- 
perature increases, Ao(T) decreases monotonically. The 
critical temperature, in terms of the chain Bethe approx- 
imation, is given as the zero of An(T). For J' / J = 0.1, 
we obtained T^/J = 0.19084(5). 



III. CRITICAL TEMPERATURE OF QID 
QUANTUM HEISENBERG MODELS 



{T,h,hcs) ^ -'^mfi,aiT,h,hcs), (14) 



where mn,a{T^h,hcs) is the staggered magnetization 
density of the a-th chain: 



mn,a{T,h,hos) = j y^(-l)''(5a,fc' 



a- 



(15) 



Here (• • • )a denotes the expectation value of the chain 
Bethe cluster and L the number of spins in the chain 
direction. 

In Fig. [21 the h^g dependence of the both sides in 
Eq. is demonstrated for J' = O.IJ and ft = 0, where 
the magnetization is calculated by means of the quantum 
Monte Carlo (QMC) method for an L = 64 chain Bethe 
cluster (see Sec. lIII Al for simulation details). It is clearly 
seen that Eq. has only one trivial solution, ftcff — 0, 
at high temperatures. On the other hand, at low tem- 
peratures two more nontrivial solutions (ftes ^ 0) appear, 
which correspond to the the symmetry-broken phase. 

The critical temperature in the framework of the chain 
Bethe approximation is defined as the point where the 
three solutions at low temperatures get degenerated with 



A. Numerical method 

In this section, we discuss the interchain coupling de- 
pendence of the critical temperature for the S — 1/2 
Heisenberg model on a QID simple cubic lattice [Eq. (|T])] 
in terms of the chain Bethe approximation. Since the 
boundary-field susceptibility of the chain Bethe cluster 
[Eq. P7)) ] can not be evaluated analytically, one needs to 
introduce some reasonable approximation or some nu- 
merical method. In the present paper, we adopt the 
continuous- time loop cluster QMC algorithm, which is 
one of the most effective methods for simulating unfrus- 
trated quantum spin systems. It is a variant of the 
world-line QMC method, based on the Suzuki- Trotter 
path-integral expansion. The continuous-time loop al- 
gorithm, however, works directly in the imaginary-time 
continuum, and thus it is completely free from the time 
discretization error. Furthermore, the correlation be- 
tween successive spin configurations on the Markov chain 
is greatly reduced, often by several orders of magnitude. 
This is manifested by the fact that clusters of spins, called 
loops, whose linear size corresponds directly to the length 
scale of relevant spin fluctuations, are flipped at once in 
the loop algorithm. 
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FIG. 3: J'-dependence of the critical temperature calcu- 
lated by chain Bethe theory (diamonds), chain Weiss theory 
(dashed lines), Tyablikov approximation (dotted lines), and 
full-3D QMC (filled squares), for antiferromagnetic (right) 
and ferromagnetic (left) interchain coupling. The QMC re- 
sults for the antiferromagnetic interchain coupling are taken 
from Ref. 0- The chain Weiss and Tyablikov approximations 
both produce the identical results for the antiferromagnetic 
and ferromagnetic cases. The solid lines are the guide to the 
eye. The error bar of each data is much smaller than the sym- 
bol size. The results of modified chain Bethe approximation 
with Jcft = 0.872J' (see the text) are also presented by open 
circles. 

The boundary-field susceptibilities (I17p are calculated 
by means of the improved estimator.^® The largest sys- 
tem size in the chain direction is L = 64 for \ J'\/ J > 0.05. 
For smaller |J'|'s, longer systems (e.g., L = 512 for 
|J'|/J = 0.01) are needed to obtain the susceptibilities 
in the thermodynamic limit L ~> oo, since the critical 
temperature decreases as |J'| does. Periodic boundary 
conditions are imposed along the chains. We interpolate 
the QMC results at various temperatures by a polyno- 
mial, and estimate the zero of Afi{T), which gives the 
critical temperature. (See e.g., the inset of Fig. 



B. Comparison of chain Bethe theory with other 
methods 

The accuracy of the present chain Bethe approxima- 
tion is checked by comparing with the results of QMC 
simulations of full-3D systems, i.e., simple cubic lattice 
of 

Lx ^ Ly X sites with periodic boundary conditions 
along all the lattice axes. For the antiferromagnetic in- 
terchain interactions, J'/ J > 0, we take the full-3D QMC 
results from Ref. 0- For J'/ J < 0, the critical temper- 
ature is estimated from the crossing point of the Binder 
cumulant 

^ ^ (19) 



TABLE I: J'-dependence of the critical temperature for anti- 
ferromagnetic (J' > 0) and ferromagnetic (J' < 0) interchain 
interactions obtained by the full-3D QMC calculation. The 
results for J' > are from the previous QMC study.— The 
figure in the parenthesis denotes the error in the last digit. 



\J'\/J 


Tc/J 






J' > 


J' < 


1 


0.94416(9) 


0.87330(2) 


0.5 


0.59248(6) 


0.55419(3) 


0.3 


0.4151(2) 




0.2 


0.30202(5) 


0.28685(2) 


0.1 


0.16917(2) 


0.16295(2) 


0.05 


0.09129(4) 


0.08900(3) 


0.02 


0.039432(7) 


0.03899(2) 


0.01 


0.020763(8) 


0.02080(2) 



of the generalized magnetization density 

for different system sizes. The largest system size we 
simulate is (L^, Ly, L,) = (40, 40, 40) and (12, 12, 264) for 
\J'\/J = 1 and 0.01, respectively. The results of finite- 
size scaling are summarized in Table ID 

In Fig. [3l we show the J'-dependence of the critical 
temperature, calculated in terms of the chain Bethe ap- 
proximation, both for antiferromagnetic (J' > 0) and 
ferromagnetic (J' < 0) interchain interactions, together 
with the full-3D QMC results. The intrachain interaction 
is antiferromagnetic (J > 0) in either case. In Fig. O we 
also present the results of the chain Weiss approxima- 
tion and the Green's function method with Tyablikov 
decoupling M'" ^ " As for the chain Weiss theory, we use 
the susceptibility of genuinely one-dimensional antiferro- 
magnetic chain calculated by the QMC method for the 
antiferromagnetic chain of 256 spins. 

In both of antiferromagnetic and ferromagnetic inter- 
chain coupling cases, the critical temperature decreases 
monotonically as \J'\/J does. For |J'| « J, it is ob- 
served that the chain Weiss theory overestimates the 
critical temperature greatly. This is is because this ap- 
proximation ignores the interchain spin fluctuations com- 
pletely. On the other hand, for | J'| ^ J the Tyablikov 
approximation gets worse. Indeed, it predicts the critical 
temperature proportional to the square root of |J'|/J, 
which does not agree with the correct asymptotic behav- 
ior, Tc ~ \ J'\/J (with some logarithmic corrections) 4*^ 
It should be noted that the chain Weiss theory and the 
Tyablikov approximation both give the identical critical 
temperature dependence regardless of the sign of the in- 
terchain coupling. 

The results of the chain Bethe approximation are fairy 
well in the whole region (Fig. [3]) . The relative errors from 
the full-3D QMC values are about 7-10% at | J'|/J = 1, 
which should be compared with the conventional chain 
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Weiss results (48-60%). Surprisingly, the present chain 
Bethe approximation predicts a lower critical tempera- 
tures for the ferromagnetic interchain coupling than the 
antiferromagnetic case, e.g., T^/J = 1.017 and 0.958 for 
J'/ J — 1 and -1, respectively. This result is seemingly 
counter-intuitive, since quantum fluctuations, which is 
expected to suppress the classical ordering, are generally 
much stronger in the antiferromagnetic cases. However, 
it is not an artifact by our approximation. Indeed, a lower 
critical temperature for the ferromagnetic interchain cou- 
pling is also confirmed by our full-3D QMC calculation 
(Table IJ). Thus, the chain Bethe approximation predicts 
not only quantitatively more accurate critical tempera- 
tures, but also its nontrivial dependence on the sign of 
the interchain coupling. 



C. Effective interchain interaction 

In order to evaluate the accuracy of the present theory 
for small | J'|/ J in a more systematic way, next we discuss 
the effective interchain coupling, J^fj(J'). This quantity, 
firstly introduced in Ref. 0, is defined as the coupling 
constant which predicts the true critical temperature if 
it is used in the self-consistent equation instead of the 
original J'. For the chain Weiss approximation, J^g is 
explicitly obtained from Eq. pT|) as 



sign(JO 
^Xc(Tc(J'))' 



(21) 



where we use the full-3D QMC results for Tc(J') (Ta- 
ble, m, which is considered to be exact within the error 
bar. On the other hand, for the chain Bethe approxi- 
mation, we obtained J^g by solving Eq. (jl6p numerically. 
In fig. m we plot the J'-dependence of the effective in- 
terchain coupling for chain Bethe approximation (J' > 
and S = 1/2) together with those of the chain Weiss ap- 
proximation (J' > and S = 1/2, 3/2, and oo)ji In both 
cases, J^ff/J' converges to a finite value for J'/ J <^ 1. 
The limiting value is 0.872 and 0.695 for the chain Bethe 
and Weiss approximations, respectively. The larger (or 
closer to unity) value of J^g/J' in the former supports 
that it is indeed a better approximation compared with 
the latter. Interestingly, the chain Bethe approximation 
has the largest J^g/ J' at J' / J = 1, in other words, 
one may say that it becomes the most accurate in the 
isotropic limit from the viewpoint of the renormalized 
factor of interchain coupling. 

On should note that, in the chain Bethe approxima- 
tion, not only the renormalization factor J^g/ J' is im- 
proved very much, but also it converges to its limiting 
value quite rapidly (already converged at J'/ J = 0.5). 
This result suggests the critical temperature might be 
well described for any values of J'/ J by the chain Bethe 
approximation using the renormalized interchain cou- 
pling constant, 0.872J', instead of the bare interchain 
coupling, J' (modified chain Bethe approximation). In 
Fig. [21 we also plot the result of the modified chain Bethe 




f /J 

FIG. 4: J'-dependence of the renormalization factor J^ti/J' 
for chain Bethe (filled diamonds) and chain Weiss (S = 1/2 
open diamonds, S — 3/2 circles, S = oo squares) approxima- 
tions. The horizontal dashed lines denote their limiting values 
(0.872 and 0.695, respectively) for J'/J 0. 



theory by open circles, which satisfactorily agrees with 
the true critical temperature in the whole range of | J'|/ J. 
Note that in this plot we use the same renormalization 
factor, 0.872, both for the antiferromagnetic and ferro- 
magnetic interchain coupling cases. 



IV. SITE-DILUTED QID HEISENBERG 
ANTIFERROMAGNET 

In the previous section, we see that the critical temper- 
ature of QID system is quantitatively improved greatly 
by the chain Bethe approximation. In this section, we 
discuss a more nontrivial example, the site-diluted QID 
Heisenberg antiferromagnet, where the chain Weiss ap- 
proximation fails even qualitatively. 

The Hamiltonian of the site-diluted Heisenberg anti- 
ferromagnet is defined as follows: 



+ J' (i,j,k^i+l,j,k^i,j,k ' S'i+ij-^fc 

~l" J ^ ^ ^i,j,k^i,j + l,kSi,j,k ' 



(22) 



where {eij,fc} are the quenched dilution factors. They 
take either 1 (occupied) or (vacant) independently with 
probability (1 — x) and x, respectively, with x {0 < x < 
1) being the concentration of vacancies, or nonmagnetic 
impurities. 

The ground state of the classical site-diluted spin 
model is equivalent to the site-percolation problemj^ 
The system undergoes a second-order phase transition 
at the percolation threshold Xp, above which there exist 
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FIG. 5: a;-dependence of the Neel temperature of the 5=1/2 
diluted Heisenberg antiferromagnet with J' j J = 0.5 obtained 
by the full-3D QMC calculation (open circles), the chain Weiss 
approximation (open diamonds), and the chain Bethe approx- 
imation (solid squares) . The percolation threshold of the sim- 
ple cubic lattice, Xp ~ 0.688, is indicated by the arrow. 

no infinite-size clusters. For the simple cubic lattice, the 
percolation threshold is determined as 

a;p = 0.6883923(4) (23) 

by the most recent simulation. ^-'^ In the quantum spin 
cases, whether a long-range order exists or not near the 
percolation threshold is a nontrivial problem due to the 
presence of quantum fluctuations. However, it has been 
established by the extensive QMC simulation that the 
staggered magnetization persists up to the percolation 
threshold on the two-dimensional square lattice. We 
expect this is also the case for the present anisotropic 
simple cubic lattice. 

In Fig.O we show the x-dependence of the critical tem- 
perature for J' I J — 0.5 (antiferromagnetic interchain in- 
teraction) obtained by the chain Weiss and chain Bethe 
approximations together with the results of the full-3D 
QMC simulation. The largest system size used in the 
full-3D QMC simulation is {L^,Ly,L,) = (16,16,64). 
The Neel temperature is estimated from the crossing 
point of the Binder cumulant of the staggered mag- 
netization for different system sizes. From the full-3D 
QMC calculations, we thus confirm that the Neel tem- 
perature remains finite at least up to a; = 0.6. For 
larger impurity concentration, it is not very easy to es- 
timate the critical temperature with satisfactory accu- 
racy in the present scale of simulation. The staggered 
susceptibility of one-dimensional chain used in the chain 
Weiss approximation is also evaluated by means of the 
QMC method. We simulate chains with spins L = 128 
and 256 for various impurity concentrations and confirm 
that there is no significant systematic difference in the 
results for those two system sizes. The Neel tempera- 
ture is then estimated by solving the self-consistent equa- 
tion (jlip numerically. Since the staggered susceptibility 
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T/J 

FIG. 6: Temperature dependence of TAn/J [Eq. HH)] of the 
S = 1/2 diluted Heisenberg antiferromagnet with J'/ J = 0.5 
a,t X — 0.40 (squares), 0.57 (circles), and 0.70 (diamonds) 
obtained by the QMC method for the five-chain cluster of 
length L — 32. The error bars are smaller than the symbol 
size. 

of one-dimensional chain diverges monotonically as the 
temperature decreases irrespective of impurity concen- 
tration, the self-consistent equation always has a solution 
(see the discussion below). 

For the chain Bethe approximation, we evaluate An (T) 
[Eq. (|16p ] by calculating the boundary-field susceptibili- 
ties of five-chain clusters with L = 32 and 64 by means of 
the QMC method. We observe any significant differences 
between the L = 32 and 64 results for < a; < 0.6 in 
the temperature range we simulate (T/J > 0.01). We 
find that the function An {T) tends to positively or nega- 
tively diverge depending on the impurity concentration. 
To see the tendency at low temperatures more clearly, we 
plot TAn(r)/J, instead of Ao(r) itself, as a function 
of temperature in Fig. [6l At low enough temperatures 
the quantity tends to converge to a finite value, which 
gives the coefficient of Curie-like behavior of Ao(T). It 
is clearly seen that for x < 0.57, the coefficient is positive, 
and thus the self-consistent equation (flB)) has a solution, 
while the coefficient is negative and Aq{T) has no zero 
for X > 0.57. 

In all the cases, the Neel temperature decreases mono- 
tonically as X increases as shown in Fig. [5] However, we 
emphasize that the result of the chain Weiss approxima- 
tion is qualitatively different from the others for large x; 
it predicts nonvanishing Neel temperature for any a; < 1, 
though the others has a finite critical concentration of 
impurities {xc — Xp and a^c — 0.57 for QMC and the 
chain Bethe approximation, respectively). Since no long- 
range order can persist for a: > a;p, the result of the chain 
Weiss theory in this regime is unphysical at all. 

Indeed, the asymptotic behavior of Tc near x = 1 can 
be discussed more precisely as follows: In the chain Weiss 
theory, the staggered susceptibility of the purely one- 
dimensional chain appears in the self-consistent equa- 
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tion (fTTj) . The percolation threshold of a single chain is 
unity, i.e., the chain is decoupled into a set of finite-length 
segments immediately by an infinitesimal impurity den- 
sity. Thus the staggered susceptibility can be expressed 
as a weighted average of contributions from finite-length 
segments: 



Xc(T) = ^p,X^(T), 



(24) 



where Xi{T) is the staggered susceptibility of a finite seg- 
ment of length £ and p£ = (1 — xYx^ the average number 
of segment of length t per site. For {1 — x) <g; 1, only 
single-site clusters — 1) contribute to the susceptibil- 
ity: 

Xc{T) = {l~x)^ + 0{{l^xf). (25) 

By solving the self-consistent equation ITTI) . the critical 
temperature is then obtained as 



{l-x)'-^ + 0{{l-xf). 



(26) 



This expression gives the exact asymptotic behavior of 
Tc of the chain Weiss theory in the vicinity of a; — 1. 
For J'/J = 0.5 and x = 0.9, Eq. ^ gives Tc/J = 0.05, 
which agrees fairly well with the result of the chain Weiss 
approximation, Tc/J = 0.043. 

A similar discussion applies also to the chain Bethe 
approximation. For x « 1, only single-site clusters con- 
tribute to the susceptibility. Since in the chain Bethe 
approximation the effective field is applied only on the 
side chains, a single-site cluster on the central chain does 
not feel the effective field, and thus the boundary-field 
susceptibility vanishes in the lowest order: 



Xnfl = 0{{l-xf). 



(27) 



On the other hand, the susceptibility of boundary spins 
is given by the same expression as in the chain Weiss 
approximation: 



Xn,c. = {l-x)^ + 0{{l~xf) fora = l,-.. ,4. 



(28) 



If these two expressions are substituted into the self- 
consistent equation (|16p . one immediately finds that it 
has no solution for < a; < 1. If one further consid- 
ers contribution from dimers (i.e., clusters consist of two 
sites) the boundary-field susceptibilities are calculated as 



Xnfl = {l-xf^ + 0{{l-xf) 



(29) 



1 



Xn,«^(l-.)^ + (l-.)^^ + 2(l-.) J 
+ 0{{l-xf) for a = I,-- - ,4. 



Again the self-consistent equation has no solution for 
J'/J > 3/4. On the other hand, for a smaller J' 
{J'/J < 3/4), there exists a solution: 



1 



2(1 



-)(! 



(31) 



However, this solution is unphysical, since it diverges as 
X ^ 1. We infer that even how one takes higher-order 
contribution from large clusters into account, there ex- 
ists no physical solution of the self-consistent equation. 
This implies that the chain Bethe approximation has a fi- 
nite critical threshold < 1, above which no long-range 
order appears at finite temperatures. 

Before closing this section, we briefly mention the ini- 
tial reduction rate of the critical temperature: 



R = 



d\ogT,ix) 



dx 



(32) 



From the present QMC results, this quantity is estimated 
as i? = 1.61 for J'/J = 0.5, which is significantly larger 
than that of the isotropic cubic lattice (i? = 1.22 and 
1.36 from renormalization-group2^ and series^ studies, 
respectively). Accordingly, the critical temperature is a 
convex function of the impurity concentration, which is 
in a sharp contrast to the linear behavior observed in 
the isotropic cubic lattice.— Such a large initial reduc- 
tion rate is also observed experimentally in quasi-two- 
dimensional Heisenberg antiferromagneti2^ The enhance- 
ment in the initial reduction rate and the convexity might 
be attributed to the spatial anisotropy of the lattice. 



V. SUMMARY 

In this paper, we proposed a novel chain Bethe theory 
for QID quantum magnets. In the present approxima- 
tion, the self-consistent equation is written in terms of 
the boundary-field magnetic susceptibilities of a multi- 
chain cluster instead of a single chain. Not only the cor- 
relations along the chains, but also those between the 
nearest-neighboring chains are taken into account ex- 
actly. As a result, the accuracy of the critical temper- 
ature of the QID Heisenberg models is improved greatly 
compared with the conventional chain Weiss theory. It is 
also demonstrated that our new approximation can pre- 
dict nontrivial dependence of critical temperature on the 
sign of interchain coupling as well as on the impurity 
concentration in randomly diluted QID Heisenberg mag- 
nets. The conventional chain Weiss approximation takes 
the random average in each chain before the thermal av- 
erage on the whole lattice, whereas the present theory 
can take fluctuations due to the randomness between the 
neighboring chains effectively. This difference in the or- 
der of thermal and random averaging has a great impact 
especially in the system with strong quenched disorder. 

In the present study, we restricted ourselves to the 
nearest-neighbor spin models on the simple cubic lattice. 
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This is because unbiased high-precision full-3D data, by 
which the accuracy of the new theory has been checked 
quantitatively, are available only for such unfrustrated 
models. It should be emphasized that, however, with 
the help of other numerical methods specialized to one- 
dimensional systems, such as the exact diagonalization 
and the density-matrix renormalization group method, 
the present chain Bethe theory itself can be applied 
straightforwardly to spin models with strong frustration 
or even to fermionic models. In such models, effects of 
correlations between neighboring chains are much more 
important, and thus the improved chain mean- field ap- 
proach formulated in the present paper could be an es- 
sential tool to investigate exotic phase transitions as well 
as anomalous low-energy properties. 
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